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The  normal  mode  solution  of  sound  transmission  has  been  used  as 


the  basis  of  an  approximate  (WKB)  calculation..  The  attenuation  of  low 

„■  .•  • -5/3'  \ 

frequency  signals  in  a surface  duct  is  shown  to  be  proportional  to  f ' > 

where  f is  the  frequency.  The  high  frequency  attenuation  is  based  upon 

experimental  data  and  is  proportional  to.  f / . vThe  leakage  of  sound 

energy  from  the  duct  is  used  as  a means  to  estimate  the  signal  received 

beneath  the  duct.  The  theoretical  results  have  compared  with  theoretical 

data.  The  surface  duct  attenuation  formula  fits  the  experimental  data 

within  the  spread  of  the  data. 
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Introduction 


The  theories  needed  to  describe  sound  propagation  in  the  half 
channel  and  surface  duct  are  essentially  the  same,  so  we  have  combined 
the  discussions.  The  ray  pechs  for  these  cases  are  shown  on  Fig.  1.  At 

very  Wgh  frequencies  when  the  surface  duct  is  hundreds  of  wavelengths 

> *■ 

deep,  the  ray  theory  gives  a good  description  of  signal  transmission  and 

* ' ; 

the  formation  of  shadow  zones.  At  low  frequencies  when  the  duct  (or 

■«  **  ’ > 

half  channel)  is  a few  wavelengths  thick,  the  wave  diffraction  effects  are 
important.  In  a shallow  duct,  low  frequency  waves  may  not  be  trapped 
and  thus  have  very  high  attenuation  while  a high  frequency  signal  would  be 
trapped  and  have  low  attenuation. 


Fig.  1 Half  channel  and  surface  duct.  The  shadow  is  crosshatched. 


A considerable  body  of  theory  and  experiment  has  shown  that  the 
depth  of  the  duct,  the  frequency  and  sound  velocity  gradients  are  needed 
to  describe  the  sound  transmission.  Since  we  must  consider  frequencies 
lower  than  have  been  studied  to  any  extent,  our  theory  should  be  extrap- 
olatable. 

The  answers  to  the  following  questions  are  needed  to  describe 
sound  transmission  in  a half  channel  or  surface  duct. 

1.  When  are  signals  trapped? 

2.  How  much  of  the  energy  is  trapped? 

3.  What  is  the  attenuation? 

4.  What  happens  to  the  energy  lost  from  the  duct  and  can  it 
be  received? 

We  have  tried  several  approaches  to  the  problem:  First,  ray 
tracing  procedures  coupled  with  estimates  of  attenuation  were  tried; 
however,  the  results  were  not  good  and  too  much  empirical  adjustment 
was  required.  Wave  theory  in  the  form  of  an  approximate  normal  mode 
calculation  was  found  to  give  agreement  with  the  experimental  data  of 
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Pederson  and  Gordon  (1965) . The  development  Is  given  in  Appendix  A.  The 
description  of  the  simplified  surface  duct  model  in  the  following  sections  is 
an  empirical  adaptation  of  those  results. 


II 


Trapped  Sound  Energy 


The  model  is  sketched  on  Fig.  2.  The  duct  is  assumed  to  be  many 
acoustic  wavelengths  deep.  More  exact  specification  can  be  obtained  from 
the  normal  mode  solution  of  the  problem  as  given  in  Appendix  A.  A slight 
rearrangement  of  equation  A20  yields  for  the  lowest  mode  the  following: 
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z < h 


where 


f 

h 

c. 
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is  the  frequency  in  cps 
is  the  mixed-layer  depth 
is  the  sound  velocity 

is  the  gradient  of  sound  velocity  in  the  mixed  layer 


On  assuming  that  part  of  the  sound  energy  is  trapped,  the  mean 
square  pressure  can  be  expressed  as  a cylindrically  spreading  wave.  With 
the  aid  of  Fig.  2,  it  is  the  following: 


= p ^ A (rh)  * A 9 exp  - Br 


(1) 


where  p^  is  the  acoustical  pressure  at  a unit  distance.  For  convenience, 
A and  B are  regarded  as  being  empirical  constants  which  can  be  determined 
as  a function  of  frequency,  h and  sound  velocity  g:adients.  Actually  an 
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Fig.  2 Surface  duct  as  a wave  guide;  x^  is  the  distance  required  for  the 
steepest  turning  ray  to  reach  depth  h.  gj  and  g£  are  the 
sound  velocity  gradients  above  and  below  depth  h.  The  effective 
depth  he  is  also  shown. 


approximate  wave  theory  solution  of  the  problem  yields  the  functions  for 
both  A and  B . 

If  the  primary  interest  is  the  surface  duct  propagation  then  the 
wave  equations  in  the  summary  of  Appendix  A should  be  used.  For  order  of 
magnitude  estimates  of  surface  duct  transmission  an  empirical  value  of  A is 
satisfactory. 


Ill 


Attenuation 


The  measured  attenuation  of  signal  transmission  in  the  surface  duct 
has  been  shown  to  have  large  attenuations  at  low  frequencies.  These 
losses  are  caused  by  poor  trapping  or  leakage  of  sound  energy  from  the 
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duct.  The  attenuation  also  increases  at  high  frequencies  and  this  loss  is 
assumed  to  be  caused  by  the  scattering  of  sound.  On  separating  the 
effects,  B can  be  expressed  as  the  following: 


B * 2 (6‘  + 6")  (2) 

where  6'  is  the  leakage  attenuation  (low  frequency) 

6"  is  the  scattering  attenuation  (high  frequency) 

Experimental  measurements  of  B have  helped  but  theoretical 
description  of  the  functional  U m of  6'  is  needed.  Diffraction  or 
leakage  of  energy  can  be  estimated,  with  the  aid  of  wave  theory. 

Normal  mode  descriptions  of  the  wave  field  were  used  by  Pederson 
and  Gordon(1965)  and  more  recently  Bucker  and  Morris  (1967) . One  could 
use  their  methods , however , we  felt  that  amount  of  computational  effort 
could  be  too  large.  A less  accurate  theory  having  simpler  equations  would 
be  sufficient  for  our  purpose.  On  the  basis  of  these  excellent  comparisons 
of  theory  and  experiment,  we  have  chosen  to  use  the  normal  mode  method 
to  develop  the  attenuation  function  B.  In  the  development  we  will  use 
their  results  as  a guide  in  making  many  of  our  approximations  and  simplifi- 
cations. The  development  and  equations  are  given  in  Appendix  A. 


The  attenuation  due  to  leakage  through  the  bottom  of  the  surface 
duct  is  (Appendix  Eq.  A2  and  A37) 


wnere 


,-5/3  , -3  -1/3  ,.6.  -1 

» f h g„  x 10  ir:  m 

6 L 


(3) 


f is  frequency  in  cps 

h is  effective  layer  depth  in  m 
e 

g is  the  sound  velocity  gradient  beneath  the  duct 

M 

h is  much  greater  than  X , and  generally  of  the  order  of  10X  . 


and 


In  a duct  having  a positive  sound  gradient,  hg  is  approximately  2h,  Fig.  2. 
Both  hg  and  the  mode  number  (here  it  is  one)  enter  es  somewhat  arbitrary 
factors  in  Eq.  (2).  Thus  it  is  reasonable  to  add  an  adjustable  factor  to 
Eq.  (2)  if  necessary. 

To  this  attenuation  we  should  add  the  losses  that  have  been 
presumed  to  be  caused  by  sound  scattering  at  the  surface.  The  high 
frequency  sound  can  also  be  scattered  by  the  smaller  inhomogenities  that 
are  generally  found  in  mixed  layers.  We  use  the  following  empirical 
formulas  given  by  Marsh  and  Shulkin  (1967): 

for  sea  state  < 3 

i /2 

a =4.5  (f/h)  ' db/kyd 

s 

for  sea  state  > 3 

l /? 

ag  = 9 (f/h)  dbAyd  (4) 

f in  kcps 
h in  ft 

or  for  exponential  function , 

6"  « (f A) 1//2  x 10  5 nepers/m 

for  f in  cps  and  h in  m (5) 

In  the  appendix,  we  estimated  the  surface  reflection  loss, 

Eq.  (A39) , and  found  it  to  be  proportional  to  frequency  squared.  The  role 
of  the  various  scattering  mechanisms  in  the  mixed  layer  are  not  understood 


and  for  the  present,  it  is  recommended  that  the  empirical  functions  (4)  or 
(5)  be  used.  The  total  attenuation  is 


B = 2 (6’  + 6") 

B - 2[f5/3h3-3g2-1/3 


x 106  + j3}  (£/h) 1/2  x 10'5] 


Use  3 for  sea  state  < 3 


6 for  sea  state  > 3 


In  normal  use  B is  converted  io  db/kyd  and  the  parameters  are  entered  in 
kH^  and  ft.  On  replacing  hg  by  2h,  the  attenuation  is 

Bdb  ■ 4 * 105  Cs)  h(£0  921/3  + {V} (h)  1/2  db/kyd  (7> 

for  sea  states  < 3 and  > 3 

The  dependence  of  B upon  frequency  is  sketched  on  Fig . 3 . The  numerical 
constants  in  Eq.  (6)  are  approximate  and  can  be  adjusted  empirically. 


Fig.  10.  Attenuation  rate  « frequency  for  the  leakage  arrival 
for  one  refraction  cycle  at  Bear,  shallow  hydrophone.  Spectrum 
levels  are  given  for  a band  width  of  1 cps  and  a reference  pressure 
level  of  l d/cm’  rms  for  1 sec.  The  height  of  the  crosses  gives 
the  variation  of  the  levels  and  the  width  of  the  crosses  gives  the 
band  width  of  the  filter.  The  theoretical  attenuation  rate  w 
frequency  was  calculated  for  a channel  depth  of  290'  and  a 
thermocUne  of  0.147  sec. 

Fig.  3 Surface  duct  attenuation,  B.  The  data  is  from  T.  Arase,  J.  Acoust 
Soc.  Am.  31,  588-595  (1959). 


IV 


Half  Channel 


The  half  channel  condition  generally  refers  to  a deep  ocean  having 
the  sound  velocity  minimum  near  or  at  the  surface . The  sound  trapping 
condition  applies  as  follows 

3/2 
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and  the  water  is  of  the  order  of  hundreds  of  wave  lengths  thick.  The 
simplest  approximation  is  to  assume  that  all  sound  energy  that  does  not 
interact  with  the  bottom  is  trapped  and  spread  cyclindrically.  Since  the 
wave  guide  is  very  thick,  leakage  losses  can  be  neglected  and  the  signal 
level  is 


.-1. 


A(rh)  A 8 exp  - ar 


where  a is  the  sound  absorption  loss  in  sea  water.  In  this  formula,  the 
vertical  distribution  of  sound  energy  is  approximated  as  being  constant 
whereas  it  actually  decreases  as  a function  of  depth.  For  more  accurate 
calculations,  the  normal  mode  solutions,  Eq.  (A8)  and  (A31),  should  be 
used. 

If  the  water  is  too  shallow  for  the  trapping  of  sound,  then  the 
problem  should  be  treated  as  a shallow  water  wave  guide.  Here  the 
properties  of  the  bottom  are  very  important. 
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Surface  Duct  Leakage 


Sound  energy  can  be  scattered  into  and  out  of  the  surface  duct. 
The  only  losses  aside  from  the  absorption  are  those  caused  by  scattering 
and  diffraction.  As  shown  on  Fig.  4,  the  flow  of  energy  is  horizontal  and 
the  losses  flow  downward. 
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Fig.  4 Surface  duct  energy  leakage 


The  energy  flow  per  unit  area  is  sound  pressure  times  the  particle 
velocity.  If  the  ray  paths  are  in  the  direction  of  energy  flow  then  the 
energy  flow  per  unit  are  proportional  to  pressure  squared.  For  our  approxi- 
mation we  will  use  the  pressure  squared. 

The  energy  loss  from  the  length  of  guide  Ax  is 

P^Ax  = [p2  (x)  - p2  {x  + Ax)]  h (8) 

2 

The  T'  ylor  expansion  for  p~  (x  + Ax)  is 

p2  (x  + A x)  = p2  (x)  + T7  P2  (x)  Ax  (9) 

O X 

and 


SbIm. 

9 x 


(10) 
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2 2 
The  loss  of  p per  unit  length  is  proportional  to  p and  is 


Afi!  --Bp2 
Ax  " Bp 


(ID 


Thus  the  leakage  pressure  is 


2 2 
Pf  = Bh  pZ 


(12) 


The  transmission  loss  in  going  from  a pressure  in  the  duct  to  a 
receiver  beneath  the  duct  is 


TL  for  leakage  - -10  log  Bh 


(13) 


The  leakage  signal  is  large  for  large  duct  attenuation  and  small  for 
low  duct  attenuation.  At  large  distances  from  the  source  the  leakage 
pressure  is  proportional  to  the  signal  in  the  surface  duct.  From  Urick  it 
appears  that  the  transmission  losses  in  going  from  within  the  duct  to 
beneath  the  duct  are  of  the  order  of  15  db.  In  the  model  that  we  have 
assumed,  the  leakage  signal  is  essentially  independent  of  the  receiver 
depth  beneath  the  duct.  The  direction  of  the  leakage  arrival  is  approxi- 
mately that  of  an  initially  horizontal  ray  just  beneath  the  duct. 

The  leakage  signals  are  traveling  in  the  forward  direction.  Thus, 
we  can  regard  the  leakage  as  being  a set  of  directional  sources  just 
beneath  the  sound  duct.  As  suggested  by  Hersey  and  Officer,  the  sound 
can  follow  the  usual  deep  transmission  paths  and  also  be  received  at 
refraction  or  convergence  zone  distances.  The  transmission  paths  are 
shown  on  Fig.  5.  Experiments  have  shown  that  the  apparent  broadening 
of  the  convergence  zone  (or  "bounce"  in  this  paper)  can  be  ascribed  to 
leakage  signals.  Fig.  6.  The  leakage  signals  following  the  convergence 
zone  were  about  20  db  less  than  the  main  signal. 
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Appendix  A 

SOUND  TRANSMISSION  IN  A HALF  CHANNEL  AND  A SURFACE  DUCT 


The  mathematical  developments  of  the  normal  mode  solution  of  the 
half  channel  and  surface  duct  are  a bit  involved  so  at  the  outset  we  will 
state  what  we  are  after.  In  a wave  guide  the  pressure  field  is 


oc 


10_ 

rh 


-6'  r 

e in 


(Al) 


We  would  like  to  calculate  the  constant  of  proportionality.  But  it  is  even 
more  important  to  derive  the  dependence  of  the  attenuation  on  frequency, 
layer  depth  etc.  Ignoring  the  loss  due  to  surface  roughness,  we  will  show 
that  the  (rms  pressure)  attenuation  is  the  following 


6'  » 
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14  m 


k„5/3h3  (3a//3 
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where 


k = — 

•'a  n . 


UL 
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h = layer  depth 

g g = sound  velocity  gradient  beneath  the  mixed  layer 
m = mode  number , 1 , 2 , 3 , ... 

The  most  optimistic  estimate  for  attenuation  is  to  assume  m = 1 . The 

-5/3  -3 

essential  behavior  is  that  6'  is  proportional  to  f h . 

m 
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The  model  that  we  used  to  estimate  the  surface  duct  attenuation  is 
that  of  Brekhovsklkh.  As  shown  on  Tig.  Al,  the  sound  velocity  is  constant 
in  the  mixed  layer  and  has  a sharp  gradient  beneath  the  layer.  The  sound 
i?  partially  reflected  by  the  change  of  the  sound  velocity  gradient.  It  may 
be  surprising,  but  comparison  of  the  attenuations  and  transmission  losses 
computed  by  means  of  Eq.  (A2)  are  nearly  the  same  as  those  given  by 
Pedersen  and  Gordon  (1965).  This  means  that,  except  for  very  high  frequencies, 
the  essential  part  of  the  trapping  in  the  surface  duct  is  the  reflection  at  the 
change  of  gradient. 


CQ  C(z) 


Fig.  Al  Surface  duct,  c (z)  = constant  in  mixed  layer 


The  acoustical  pressure  due  to  a continuous  wave  point  source  in 
a wave  guide  has  been  discussed  in  much  detail  by  Tolstoy  and  Clay  (1966) 
We  will  use  their  results  and  generally  follow  their  notations.  In  the 
stratified  wave  guide,  c (z) , the  acoustical  pressure  at  large  range  is 
(with  the  aid  of  their  Eq  2.5  and  3.118)  the  following: 
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= -M2-,  Y P exp  - i (x  r - wt  - ~ - 6 r 

1/2  L m m 4 m 


m = 1 


P = p Q (zj  <p  (z)  p 
m m m 0 m s 
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(A3) 


(A4) 


_ 0<T  , ni1/2  -1  -1/2 

p = 2n  (p  c-JL)  p v x 

m s 0 s m m 


(A5) 


Zq  and  z are  the  source  and  receiver  depths  (Fig.  Al) 
r is  range  in  cylindrical  coordinates 

z is  depth,  positive  downward 

U>  is  angular  frequency 

x is  horizontal  component  of  wave  number 
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m 
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k 

ji 


= a’/c  , c is  phase  velocity 
in  m 3 

~ vertical  component  of  wave  number 

= (y  2 + x 2)  1/2  = u»/c  = 2ttA 
m m 

is  source  power 
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p and  p are  the  densities  at  source  and  receiver 
s 

6 i&  the  mode  attenuation  and  is  an  addition  to  the 

m 

original  equation 

are  the  eigen  functions  and  are  solutions  of  the 
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separated  wave  equation: 


9"  A 2 . 
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where  c = cQ  (1  - 28^) 
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The  solutions  of  Eq.  (A6)  in  a stratified  wave  guide  (7  is  a function  of  z) 

m 

generally  require  extensive  numerical  computations . Our  purpose  in  the 
following  is  to  develop  approximate  solutions  which  can  be  used  as  a guide 
in  devising  a simple  approximation  such  as  Eq.  (Al) . But  first  we  take  the 
absolute  square  of  Eq.  (A3)  and  examine  the  result. 
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P p e'(6m+6n>r 
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cos  (x  - x ) r 
m n 


(A8) 


The  first  summation  gives  the  average  level  as  a function  of  range.  The 

second  set  of  terms  is  oscillatory  and  these  terms  describe  the  maximums 

and  minimums  often  observed  in  experiments . These  terms  depend  upon 

(x  - x ) and  great  accuracy  in  the  computation  of  the  eigen  functions , x , 
m n m 

is  required  to  theoretically  duplicate  the  actual  maximums  and  minimums  of 

experimental  data.  At  very  large  range,  the  sound  field  is  analyzed  by 

measuring  the  mean  level  and  correspondingly,  we  ignore  the  second 

summation  or  the  fluctuations.  Thus,  approximate  eigen  values  x will 

m 

be  adequate. 


-4- 


We  use  the  WKB  approximation  (Tolstoy  and  Clay,  pp.  48-52)  to 
determine  the  eigen  functions.  As  shown  on  Fig.  A2,  the  following  sound 
velocity  profiles  are  convenient  to  describe  the  surface  duct  and  for 
computation: 

-1/2 

c * Cq  (1  — 2a^z)  0 < z < h 

- 1 /2 

c * c,  [ 1 + 2a_  (z  - h)  ] z>  h 
n L — 

ch  = c0  (1  - 2ajh)“l/2  (A9) 

To  begin  with,  we  are  assuming  a positive  gradient  in  the  mixed  layer  so 
that  we  can  also  apply  the  theory  to  the  half  channel. 

For  small  gradients,  Eq.  (A9)  are  the  approximate  linear  functions 
as  follows: 

c » Cq  (1  + a^) 

c « ch  t 1 ~ a2  (2  “ h)  ] (A10) 

The  approximate  eigen  functions  are 

-i/2 

<p  (z)  y sin  [s  (z)  + s.]  (All) 

m o c m mi 

s (z)  = V 2 y dz  + sA  = mir 
m J m 0 

0 

m = 1,  2,  3,  . . . M (A  12) 
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Fig.  A2  Surface  duct 


The  condition  for  validity  of  the  WKB  method  is  (Tolstoy  and  Clay,  2.213) 


- -j-  2r,  y <<1 

y dz 


(A  13) 


The  approximate  eigen  functions,  Eq.  (All),  cannot  be  used  at  small  y , or 
at  the  turmng  depths.  The  mathematical  techniques  used  to  connect  the 


WKB  eigen  functions  through  turning  points  are  discussed  in  many  quantum 


mechanics  texts.  We  will  not  use  p near  the  turning  depths  and  not 


calculate  the  connection  functions. 


Modes  are  trapped  for  those  frequencies , or  x , for  which 


Eq.  (A12)  is  satisfied.  Let  us  assume  that  the  mth  mode  turns  at  depth 


z , has  the  wave  number  component  x and  phase  velocity  c . y is 
m m mm 


2 ,2  2 2 ,-2  -2. 
y = k - x = co  (c  -c  ) 
mm  m 


(A  14) 


a * ; 
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With  the  aid  of  the  preceding  equation,  Lhe  integration  of  Eq.  (A12)  [for 
c given  by  Eq.  (A9)  ] yields  after  manipulation 


m 


(z  ) = t*~  (2a. z )3//2  + 7 = mrr 

m 3a  ^ 1 m 4 


(A  15) 


where  (A9)  or  (A10)  have  been  used  to  express  the  parameters  as  follows: 
K = <o/cn 


s (z)  = 
m 


m 


(m  - “)  it  [ 1 - (1  - 2/2  )3//2  ] 
4 m 

k0  (2alzm)1/2  (1  - z/\)1/Z 


m 


k0  <1'2alzra) 


1/2 


hQ  (1  - 28J  2) 


1/2 


m 


[ (m  - ~)  ff  k0-1]2/3  [2c*l  ]_1 


2a 


(Note  in  Eq.  (A14),  lot  y =0  to  determine  c ) 
' m m 


In  summary: 

WKB  is  valid  for  (by  application  of  Eq.  (A13)  ) 


2k0zm<2alzJ1/2  (1-z/zm)3/2  » 1 
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(A16) 


(A  17) 


For  trapped  modes 


A 


fg 


(23^) 3/2  > (M  - ~ ) ir  (A  18) 

M is  an  interger 

For  a given  h,  and  a^  , k^  must  satisfy  (A18)  for  the  waves  to  be  trapped. 
The  WKB  eigen  function 

<P  (z)  oc  (1  - ~ ) 1//2  sin  Cs  (z)  - ^ 3 (A  19) 

m zm  m 4 

(Note  that  s.  = -tt/4  because  (p  vanishes  at  the  free  surface  or  z = 0) 

1 m 

The  first  question:  "When  are  signals  trapped?"  can  be  answered 
by  application  of  condition  (A18) . A rearrangement  gives 

kQh  > 3(M-^)ir  (2a1h)"l/2/2  (A20) 

An  evaluation  of  (A19)  for  the  surface  duct  of  Pedersen  and  Gordon  (1965) 
follows: 


sound  velocity  gradient  « 10  sec 


h « 100  m 


and 


aj  « 7 x 10 
M » 1 


-1 

m 


(A21) 
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Thus 


30 


koh  > 


£ = 0)/2ff  > 225  cps 


To  trap  the  first  mode,  this  duct  must  be  more  than  15 X thick.  A similar 
calculation  for  second  mode  yields  f > 530  cps.  For  530  cps  signals,  they 
show  fair  trapping  of  the  second  mode  and  large  leakage  for  higher  modes . 

Recalling  Eq.  (A3)  and  (A4) , it  is  evident  that  P is  constant  for  a 

m 

fixed  source  and  receiver  depths,  power  and  wave  guide  parameters.  P 

m 

depends  upon  the  product  <p  (zft)  <p  (z)  and  thus  upon  the  source  and 

mum 

receiver  depths.  Since  c o tends  to  zero  as  z tends  to  zero,  the  source  and 

m 

receiver  depth  dependence  should  be  included.  Incidentally,  the  value  of 
<Pm(z)  for  small  z is  the  same  as  that  given  by  the  dipcle  formula.  To  show 
this,  we  let  Tm  be  defined  with  the  aid  of  the  grazing  angle  0 as  follows: 


y = k sin  0 
m u 


(A22) 


The  substitution  of  Eq . (A22)  in  (A3)  and  (A7)  yields 


2 1 2 
p oc  sin  7 z 
' m 


(A23) 


The  evaluation  of  s (z)  for  small  z gives 
m 


<Pm  (z)  «sc  sin  y z 
m m 


for  z/z  <<  1 
m 


P ©c  sin 
m 


7 Z 

m 


(A24) 
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"Ha  


We  observe  that  the  dipole  effect  is  included  in  the  mode  formulation  and 

has  the  same  form  as  it  would  have  in  a ray  trace  solution.  Since  max  7 

1 m 

or  sin  6 can  be  determined  very  easily  for  a simple  surface  duct,  it 
max 

is  easy  to  include  the  surface  reflection  for  very  shallow  sources.  The 
comments  concerning  deep  dipoles  and  uncertain  source  depths  would  also 
apply  to  the  eigen  functions  for  deep  source  and  receiver. 

Eigen  functions  as  a function  of  depth  are  sketched  in  Fig.  A3.  The 
approximate  functions  and  exact  ones  from  Bucker  and  Morris  are  compared. 
The  limit  of  validity  of  the  WKB  approximation  is  indicated. 


Exact  [<p| 


WKB  Jp| 


Limit 


Fig.  A3  Eigen  functions 


a)  From  Bucker  and  Morris  (1967) 

b)  WKB 


Approximate  evaluation  of  the  excitation  functions  can  be  made  by 

using  the  WKB  eigen  functions  (All)  and  ignoring  the  dependence  of  7m 

upon  z.  (Vm  is  slowly  varying  except  near  the  turning  depth.)  In  this, 

the  infinite  value  of  <p  (z)  at  the  turning  depth  is  eliminated.  We  also 

m 
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approximate  (p  (z)  by  letting  tp  (z)  be  zero  for  depths  larger  than  zm . 

mm  m 

The  approximate  (pm  and  derivatives  are 


-l/2 

tp  (z)  « y sin  is  (z)  + s.  ] 0 < z < z 

m m m 1 ~ ~ m 


(A25) 


ym  » constant 


O (z)  = 0 

m 


z > z 


m 


(A26) 


3z 


« y cos  [s  (z)  + s ] 
m 1 


m 


by  Eg.  (A12) 


d s_  (z) 


m 


d z 


= y 


m 


In  Eq.  (A7),  v becomes  after  substitution  of  (A25)  and  (A26)  the  following 


I'm  * c02  (2to2a  )-1  [ym  ] 


z = 0 


(A27) 


We  can  express  the  source  function  as  either  a source  power  J”L 
n rms  pressure 
Clay  (2.46)  - (2.50). 


or  as  an  rms  pressure  at  unit  distance  R^.  The  relation  is  from  Tolstoy  ana 


2 1 ,-Tlc  , 1 

p I = ? 


tone) 1/2  = 2 d 1/2  Pq  rq 


(A28) 

(A29) 
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The  substitution  of  (A27) , (A29)  and  (A16)  into  (A5)  yields 


4 n 3/2  kQ1/2  (2ax) 

m (2a  z )1/2  (1  - 2a  z ) 1/2  P°  0 
1 m 1 m 


(A30) 


D d ,.3/2  -1  -i  . -l/2  0 .-1/2 

pn,  = p zm  kn  (1-2alzm> 


sin  [s  (z)  - “ ] sin  [ s (zj  - T ] (A31) 

m 4 m 0 4 

Numerical  values  of  P can  be  used  in  (A3)  and  (A8)  for  the  calculation  of 

m 

the  acoustical  pressure. 

Before  we  estimate  the  attenuation,  let  us  consider  our  results  in 
the  context  of  Eq.  (Al)  and  Fig.  2.  Radiation  into  the  M modes  corres- 
ponds to  signals  leaving  the  source  at  a succession  of  angles  0j,  6^, 

. . . . 0^  is  approximately  A0  in  Fig.  2.  The  factor  r-3  Pm^ 

includes  (rh)-1.  As  mentioned  earlier,  the  dipole  or  surface  reflection  is 
contained  in  sin  y [ ] . 

The  half  channel  can  be  evaluated  by  means  of  Eq.  (A8)  and  (A31). 
The  attenuation  would  be  small  and  the  measured  attenuation  of  signals 
in  deep  water  can  be  used. 

The  attenuation  in  the  surface  duct  is  caused  by  leakage  of  waves 
through  the  bottom  of  the  duct.  The  amount  of  leakage  depends  upon 
gradient  of  sound  velocity  beneath  the  duct  in  addition  to  the  parameters 
that  describe  the  duct  and  the  signal . The  velocity  gradient  below  the 
duct  is  of  the  order  of  40  times  the  gradient  in  the  duct.  For  a rough 
approximation  to  the  form  of  the  attenuation  function,  we  will  follow 
Brekhovskikh  (English  ed  1960)  and  let  the  gradient  in  the  duct  be  zero. 
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We  give  here,  simplified  results  that  are  valid  for  trapped  signals. 
The  actual  computation  is  for  the  modes  that  are  reflected  by  the  change  of 
the  velocity  gradient.  After  considerable  manipulation  between  pp  528  and 
539,  he  gives  the  following.  (Where  we  have  used  q for  his  s,  m for  his  1 


and  6 for  his  0 on  pp  537-539): 
m m 


3a_ \ 1/3 


11  5 » koh(^f) 


(A32) 


for  q >>  1 


(trapped  modes) 


m ff 

Xm  3q  + i\/3~  (1.57)/2 


(A3  3) 


x = x‘  - ’ x" 
m m m 


(A34) 


6’_  * (9a  2kft)1/3x’  x" 

L U mm 


m 


6’  * 14  m2  [ k 5/3  h3  (3a  ) 1/3  ] _1 

m u c 


(A3  5) 

(A3  6) 


Brekhovskikh  gives  a table  of  x'  and  x"  as  function  of  q (or  his  s)  for 
high  and  low  frequency,  p.  540,  and  it  is  our  Table  3.1. 

The  attenuation  of  a 530  cps  signal  in  a 100m  duct  above  a negative 
gradient  (a2  = 3 x 10~^  m-*)  is  about  4 x 10”^  m“*  or  about  0.3  db/km. 
These  attenuations  are  somewhat  higher  than  reported  by  Pedersen  and 
Gordon  but  it  should  be  emphasized  that  we  have  assumed  no  trapping  of 
energy  in  the  surface  duct  other  than  that  reflected  by  the  change  of 
gradient  of  c at  z = h. 
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A simple  way  to  approximate  the  effect  of  a duct  having  a positive 

gradient  is  to  replace  h in  Eg.  (A36)  by  an  effective  h . To  estimate  the 

effective  depth,  let  us  assume  that  the  ray  leaving  the  surface  at  0 

travels  straight,  reflects  at  the  depth  h and  returns  to  the  surface  at  the 

same  distance  as  the  actual  curved  ray.  For  very  small  0 , the  effective 

m 

depth  is  approximately  2h.  Thus  (A36)  can  be  written  as 

6‘  » 14  m2  Ck  5/3  h 3 (3a, ) 1/3  3 ~1  (A3 7) 

m 0 e 2 

The  value  of  h should  be  regarded  as  adjustable. 

We  need  to  include  a second  form  of  attenuation,  the  loss  at  the 
rough  sea  surface.  On  the  basis  of  reflection  of  signals  at  a rough  surface, 
Tolstoy  and  Clay  have  estimated  the  attenuation  in  a wave  guide  having 
irregular  boundaries.  Their  equation  (6.120)  is 

6"  = y3cr2  (h  x )-1  (A3  8) 

m m e m 

and  with  the  aid  of  (A  16) 

6 -•  kn  ct  (2a.  z ) h U-2a.  z ) 

m 0 1 nr  e 1 m 

In  approximation,  we  can  let  z and  h both  be  equal  to  h,  so  that  6"  is 

me  m 

6"  « k2  a2  (2a,h)3/2  h"1  (A3 9) 

iv-  1 

The  total  attenuation  number  is 

6 =6'  + 6"  (A40) 

m m m 
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Summary: 


For  waves  to  bo  trapped 

k0h  > f (M  " 4 (2a  j h)"1//2 

(A  19) 

| 1 « M 

p 2 . £-  V p 2 -2  6 r 

11  r L Pm  0 + cross  terms 

m =»  1 

(A8) 

, 3/2 

P ft>  f,,  o \ 4 ff 

0 o'  ~l7i  ift  sill  Ck  <2a  K ) 

k0  v(1-2ai2,„)  0 1 “ 

1/2  , 

~0J 

sin  [k  (2a. z )i//2  x ] 
u 1 m J 

(A3 1) 

- ■ 1 10ho(3a2A0)1/3 

(A  3 2) 

- Mm2  CkQS/3  ho3  (302,l/3  j -1 

(A3  6) 

i m 0 1 

in  which 

(A3  9) 

0 ■ °0  0 - 2a,  s)’1/2  0 < e < h 

o - chn  + 2„2  k-io  r1/2  a j.  h 

(A9) 

zm  “ [3ai  k0'*]2/3  [a^r1 
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a m rms  sea  surface  roughness 

h * offoctivo  depth  of  corresponding  constant  velocity 

e 

surface  duct,  hQ  » 2n. 
p is  the  source  pressure  at  distance  RQ. 


For  one  mode 


, the  transmission  loss  in  db  Is  relative  to  pQ  at  RQ 


the  following: 


log  (|p|  2/pr>2)  a 10  i0Q  (l/V 


1.2xl02.»2CKnS/3l>  3 (3a2)l/3r'r 


0 e 


+ a . 6 k02  a 2 (2a y £he  (l  ~ 2d l *n?  1 


Ar  + 10  log  (k0  «m“A0) 


j/2  1 

20  log  jsln  k()  (2a  j 2o> 


20  log  {sin  kg  (2a  j 2m) 


1/2  4 


+ 10  log  (1  - 2a1  zj  - 10  log  16" 


(A41) 


ond  whoro  A la  the  absorption  loss  lor  sound  propagation  In  sea  water. 

As  a numerical  example,  wo  assume  conditions  that  are  approximately 
those  of  Pederson  and  Gordon: 


m- 


m 


s- 

JL 

1 


1. 


6 


v; .. 


5 * 

SE  ! 


§|^#§l|f|§P?g<S 


f = 530  cps,  h = 100m, 


7 in“6  -1 

a^  - 7 x 10  m 


-4  -1 

a„  = 3 x 10  m , a = 0 , z = z = 16m 

L U 


and  ignore  A.  The  results  of  calculation  are  given  on  Table  1.  Two  values  of 
the  effective  depth  are  given  to  demonstrate  the  effect. 


Table  1 Estimates  of  mode  attenuation  61 


m 


M = 2 , kg  ~ 2m-^ , (2a^  zm)  « 4 x 10“^ , 


= 100),  q«5 
m 


Pedersen 


Mode 

Effective  depth 

Estimated  attenuation 

and  Gordon 

for  m = 1 

h = 100  m 

8.6  61  « 
m 

0.3  db/km 

h = 200m 
e 

8.6  6'  « 
m 

0.04  db/km 

7 x 10" 3 db/kyd 

for  m = 2 

h = 100m 
e 

8.6  6'  » 
m 

1 . 2 db/km 

h = 200m 
e 

8.65'  « 

m 

0.16  db/km 

0 . 3 db/kyd 

roughness 
a - 0 


6"  = 0 
m 


0 = 2m 


h = 100m 
e 


8.6  6"  « 0.2  db/km 
m 


% 

M 


Pedersen  and  Gordon  show  data  for  a relatively  short  range  of  10  kyds  and 
it  is  difficult  to  measure  the  attenuation  coefficients  of  the  low  modes  in 
this  distance.  The  estimated  attenuation  are  probably  with  the  right  order 
of  magnitude. 


Wii  f'iii  >'i.ri  "im  ^ 


?W«5Sl . _ 

- — = -JU*4 


Evaluation  of  the  rest  of  the  constants  in  Eq.  (A41)  yields  'or  m = 1, 
h^  = 100m,  a = 0,  and  r in  km,  i.e.  = 1 km 


-10  log  (|p|2/P02)  « 52  + 10  log  (r/Rj)  •*  0.3r,  ^ = 1 km  (A42) 

If  more  than  one  mode  were  contributing,  Eq.  (A8)  should  be  used.  A 
comparison  of  Eq.  (A42)  and  their  data  are  shown  on  Fig.  A4. 


iv 

4Cll 

_ ft* 


Eq.  (A42) 

" — - — _ 1 mode 


— .t.v' 


RECEIVER  DEP1H  50  FT 
FREQUENCY  5JO  CPS 

MODE  THEORY  R S PROFILE 

MODE  THEORY  S S PROFILE 
EXPERIMENT 


i -4 4 4- »— 

RANGE  IK  YDS) 

Fic.  5.  Propagation  losses  (or  50  ft  receiver,  55-tt  source. 


Fig.  A4  From  Pedersen  and  Gordon,  1965. 


In  conclusion,  we  suggest  the  next  steps.  At  high  frequencies,  the 

number  of  trapped  modes  is  high  and  correspondingly,  Eq.  (A41)  should  not 

be  used.  Brekhovskikh  has  studied  this  and  has  suggested  the  use  of 

, , 2 

"average  decay  laws",  pp.  415-421  (1960).  He  averaged  p over  the  depth 

and  then  replaced  the  summation  of  modes  by  an  integration.  For  mode 

2 2 
attenuations  proportional  to  m [Eq.(A36)]  , the  average  decay  of  p is 


proportional  to  r 


-3/2. 
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